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Abstract 

A neutral Dirac fermion ifi is supplied as a singlet within the context of the standard model (SM) 
and is considered as a dark matter (DM) candidate near electroweak scale (10 ~ 1000 GeV) with 
nonzero magnetic dipole moment. The Dirac particles have four different types of electromagnetic 
couplings (four form factors) in general. We expect that the candidate mainly interacts with SM 
particles through magnetic dipole moment (MDM), since MDM conserves the discrete symmetries 
like parity (P), time reversal (T), and charge conjugation (C) or its combination CP. The magnetic 
dipole moment constrained by the relic density may be as large as 10" 18 ~ 10~ 17 e cm. We 
show that the elastic scattering is due to a spin-spin interaction for the direct detection, and the 
predictions are under experimental exclusion limits of the current direct detectors, XENON10 and 
CDMS II. We also consider the possibility of WIMP detection in near future. 
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I. INTRODUCTION 



The identity of the dark matter (DM) has been one of the most important outstanding 
questions in physics and cosmology. It is known that DM must be nonbaryonic, nonrela- 

n 

tivistic (cold) and nonluminous matter So it is not one of elementary particles contained 
in the standard moldel (SM). The leading candidate from particle theory is a weakly inter- 
acting massive particle (WIMP), since the annihilation cross section can lead to the correct 
thermal residual abundance with the right magnitude. The new stable massive particles 
could be generated in many extensions of the SM 2] {3] {4] jsljs]. 

In this letter we consider a Dirac fermionic dark matter candidate with mass near elec- 
troweak scale (10~1000 GeV), which is supplied as a singlet within the standard model 
context. The Dirac particles have four different types of electromagnetic couplings (four 
form factors) in general. Among the four couplings, we expect that magnetic dipole mo- 
ment (MDM) is sizable for the neutral fermion. The electric charges must be prohibited 
by the unbroken £7(1) q gauge symmetry [7] and the other couplings except for magnetic 
dipole accompany violations of the discrete symmetries. The violations of the discrete sym- 
metries might be sizable for the DM, but the violations have been estimated to be small for 
the known Dirac particles. We assume that the candidate has the small violations of the 
discrete symmetries. 

The singlet Dirac fermion ip should couple to the shifted photon (hypercharge gauge 
boson B) near electroweak scale (10~1000 GeV), instead of photon. The Lagrangian with 
the magnetic dipole moment of B, is 

C e ff = i/i£^V^ (1) 

In the standard model context, the Lagrangian may be expressed about photon and Z 
boson. 

£*ff = \v4(J^(F^ - taaOwZ^) (2) 

where F^ v is the field strength for photon, Z^ u for Z boson and fi is the magnetic dipole 
moment 1 . 



1 To be strict, this is a magnetic dipole form factor, i 7 ^? 2 )- The magnetic dipole moment is defined at 
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FIG. 1: Feynman diagrams for WIMP pair annihilation into f f,WW, ZH. The hatched circle 
indicates the vertex for the dipole coupling. 

II. ANNIHILATION PROCESS 

In the leading order of the magnetic dipole moment, WIMP annihilates into SM particles 
(//> WW, ZH) via the s-channel 7, Z exchanges, where the Feynman diagrams are given 
in Fig.l. The MDM operator is split into energy and momentum dependent parts by the 
familiar Gordon identity. 

v(p')a^q v u(p) = iv(p) (2M^ + (p - pY)u(p) (3) 

The MDM operator has energy dependence. This implies that there is no threshold or 
helicity suppression due to violations of discrete symmetries. There is also no lepton number 
violation for this annihilation process, since the candidate is a Dirac particle. The below 
calculation for annihilation rates confirm this fact. 

The squared scattering amplitude in an average over spin of initial states and a sum over 
colors and spins of final states is 

\M\ 2 = |Q 7 P 7 - Q z t&n6 w P z \ 2 L^OTm™ (4) 

where Q y , Qz are couplings involved the outgoing particles with photon and Z boson, P 7 = 
1/s, Pz = 1/ [(s — m 2 z ) + im z Tz} imply photon and Z boson propagations, and N c is 
number of colors for the final states, 3 for quarks and 1 for the other particles. L^ v is the 
trace of initial states and Q pa stands for the trace of final states. 

q 2 = 0. However we do not distinguish the difference of both terminologies in this letter. 
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FIG. 2: (color online) Annihilation fraction of particles. The fermionic annihilation products 
dominate, hadronic annihilation products in particular. 



The consistant annihilation rates for each final state are 

IV/-'y 

e/P 7 



f 



N c ap f 



2 cos 2 6 



w 



(s + 2m)) + 



c A Pz 



2 cos 2 6 W 



[s-Am)) 



/x 2 (s + 8M 2 ) 



&W v re\ 



VZHVrel 



a(3 w s 
24 

a \p\ 



2 s 2 + 20sm 2 w + Y2m% / 



rn 



/i 2 (s + 8M 2 



w 



Tflr 



mw cos 2 6 



-P? 



1 + 



3m| 



•/i 2 (s + 8M 2 



(5a) 
(5b) 

(5c) 



where t> re i is the relative velocity of the annihilating particles, a is the electric fine structure 
constant, (3 is the velocity of the particle, cy = T| — 2e/ sin 2 9w, ca = Tj are the vector and 
axial couplings of Z-boson with electric charge Cf and the third component of isospin T|, 
and |p| = 277j(( s — m H + m z) 2 — 4sm|) 1//2 is momentum of Z-boson or Higgs. The total 
annihilation rate is av vc \ — (erf + (?w + ctzh) v xe \. We show the explicit mass dependence of 
fermions. 

Fig. 2 shows the annihilation fraction as a function of WIMP mass for the final states. 
The fermionic annihilation products dominate, particularly hadronic annihilation for entire 
parameter space. ZH channel has the small contribution since it interacts with WIMP via 
only Z-boson. The 7 — Z interference plays very important role for the VF-boson channel 2 , 



2 If we only consider the interaction with photon, the contribution by VV*-boson becomes extremely large in 
high energies. The interaction with Z -boson avoids this danger. 
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since the VT-boson channel contribution results in a significant decrease in high energy due 
to the interference. 



III. CONSTRAINT FROM RELIC ABUNDANCE 

In the very early universe, these particles would be present in large numbers in thermal 
equilibrium. They could reduce their density only through pair annihilation as universe is 
cool. This pair annihilation would be efficient enough to reduce the present day number 
density of WIMP. The time evolution of the relic density is described by the Boltzmann 







After the freeze out, the 



equation and we can solve the equation in the approxmate way 
actual number of the dark matter per comoving volume becomes constant and the present 
relic density is determined. The relic density is roughly 



, (1.07 x 10 9 )x F 



VtcBMh ~ — — — , , r (6) 

wg^mpi{GeV){av Tel ) 



where g* is the effective degrees of freedom in equilibrium, mp\ is Planck mass, xp is the 
inverse freeze out temperature and (av Te \) is the thermal average of annihilation rate. 

The inverse freeze out temperature has almost no dependence on WIMP mass, xp ~ 20 
near electroweak scale. The present relic density of CDM measured from WMAP (Wilkinson 
microwave anisotropy probe) data [9] is ^cdm^ 2 = 0.111±0.012 at the level of 10% accuracy. 
So the thermal average of annihilation rate is (<Jv rc \) ~ 0.62 pb. 



We take the thermal average 
tics, which is described in Ref. 



or the annihilation obtained with Maxwell-Botzmann statis- 



3. 

where Ki is a modified Bessel function of order %. 

Combination of Eq. (5) and (7) gives the allowed parameter sets for the magnetic dipole 
moment of WIMP. Fig. 3 shows the allowed parameter sets which satisfy the WMAP mea- 
surements for relic density, in the case that the Dirac fermions account for all the observed 
dark matter. So those become the upper limits of WIMP magnetic dipole moment. All the 
annihilation rates of Eq.(5) are included and the known inputs are used as well as Higgs 
mass of 120 GeV. The anomaly indicates the resonant region of the Z boson exchange where 
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FIG. 3: Allowed parameter sets from the relic abundance. The abundance is chosen to be 
ficDM^ 2 = 0.111. The arrows indicate the points where new channels start. 

2M ~ mz- In this region, the dipoles are relatively small in order to compensate the en- 
hancement of annihilation rate by Z boson resonace effect. The WIMP magnetic dipole 
mmoment may be as large as 1CT 18 ~ 10~ 17 e cm. If the dipole comes from new physics 
effect, we predict the new physics with 0(A) ~ 10 TeV. 

IV. DIRECT DETECTION 

The sensitivity of dark matter is controlled by their elastic scattering with a nucleus 
in a detector. The Feynman diagram which give rise to the WIMP-nucleon interaction is 
shown in Fig.4, which the t-channel exchanges of 7, Z are only possible. For the cold dark 
matter, the momentum transfer is quite low, \q\ ~ 0(MeV)(En ~ O(KeV)). The interaction 
via Z boson is l/m| ~ Gf where Gf is the Fermi constant, and 1/q 2 for the photon. The 
interaction via Z boson may be neglected. We thus only consider the interaction via photon. 

A. Elastic-scattering cross section 

The effective Lagrangian may be described by four fermion interaction, 

Cq = i^L (V^V) (*VZ) (8) 
where e is the electric coupling and e q is an electric charge for the quark q. 
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FIG. 4: Feynman diagram relevant to WIMP-nucleon elastic scattering. The hatched circle indi- 
cates the vertex for the dipole coupling. 

The MDM operator ipa^ q v ip has no energy dependence for the elastic scattering. So we 
take only space component of the operator and the leading order of \q\ in the nonrelativistic 
expansion of the spinor. It results in the spin-spin insteraction between WIMP and quark. 

The effective Lagrangian for nucleons may be expressed by the spin distrbution of quarks. 

C N D 2iM^- Prty) (N%N) (9) 

where the symbol iV stands for proton (p) or neutron (n), and = ^e q A^ N 'q where 
q is the fraction of nucleon spin carried by quarks of type q. The fractions for proton, 



according to Ref.pjj, are A^u = 0.78, A^d = -0.48, and A^s = -0.15, each with an 
uncertainty of about ±0.02, and A^u = A^d, A^d = A^u, A&h = A^s for neutron. 

The differential cross section of WIMP-nucleus in a spin average of initial states and a 
spin sum of final states is given by 

^e L = J^_(J_V 1 E lM] 2 (1Q) 

dq 2 47rv 2 q*\2m) {2S + 1)(2 J + l)^ 1 1 1 ' 

where S = § is the WIMP spin, J is the nuclear spin and m is mass of the atomic nucleus. 

The squared scattering amplitude can split into two parts, WIMP and nuclear parts. The 
WIMP part is 



(2S + 1)^' IW1MF 3 
The squared amplitude for the nuclear part depends on the nuclear structure function, 
S'(lgl), which is normalized by the nuclear spin. 
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^j^E\M\^ cus = ^p±S(\q\) (q 2 5- - ,V) (11) 

The nuclear structure function 5(|q|) may be expressed in terms of three form factors, 
isoscalar (do), isovector (ai) and the inteference (aoai). 

S(\q\) = a 2 Soo (\q\) + a oai S 01 (\q\) + a\S u (\q\) (12) 

with a = a p + a n and a\ = a p — a n . The functions encompass the magnitude of the spin 
associated with the nucleon populations, as well as the effects of the spatial distribution of 
that spin at nonzero momentum transfers. 

The WIMP-nucleus differential cross section results in 

where /i\ = (^J) 2 S(0) is the magnetic dipole moment of the atomic nucleus and F 2 (\q\) = 
S(\q 2 \)/S(0) is the normalized form factor, F 2 (0) = 1. We have to integrate over q 2 for the 
actual cross section. There are two quantities which depend on atomic nuclei, nuclear spin 
J and the nuclear structure function S(\q 2 \) ~ fi 2 A F 2 (\q\). Those become the reasons we 
have the different signatures for the various targets. 
The cross section at \q\ = is 

^ 2 da el (\q\ = 0) 2 2 2 J+ 1 
^ ^ = 3^ M ^ (M) 

where M r = Mm/ (M + m) is the reduced WIMP-nucleus mass. 
The WIMP-nucleon cross sections result in 

1 2 2 /ea Pi „\ 2 J + 1 
= 2^ m ^ l^-J — (15) 

where m r = Mm^l{M + m N ) is the reduced WIMP-nucleon mass. This cross section 
still includes one of the target properties, nuclear spin J, so it is not totally the nucleus 
independent quantity. 
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B. WIMP detection 



The expected event rate depends on WIMP-nucleus cross section and the WIMP flux on 
the Earth, so the differential event rate per unit target mass and unit time is 



dR 



Pd 



V f{v)—r^rdv 



(16) 



dq 2 mM J dq 2 
where pr> ~ 0.3 GeV/cm 3 is the local DM density in the solar vicinity and /(f) is the distri- 
bution of speeds relative to the dectector. We take the Maxwellian distribution considered 
the annual modulation |l2| . 
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(17) 



where vq = 220 km/s is the circular speed of the Sun around the Galactic center and 
ve ~ Vq is the Earth speed to the Sun. ve changes as the Earth's motion come into and out 
of alignment with the Sun's motion around the Galaxy. This causes an yearly modulation 
in the event rate which peaks around June 2nd each year [l3j. We take the average velocity, 
ve = 232 km/s. 

With the relation q 2 = 2mE R , Eq. (13), (16) and (17) lead to the differential event rate. 
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The final formula for the event rate per unit target mass and unit time is given by 
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(19) 



where v m w 



ER 2M^ m • T ne f° rm factor has to be integrated over E R in the recoil energy 



range. The form factors for the spin-spin interaction are relatively slowly decreased, so it 
results in suppressions for the event rates. 

There are various isotope targets which have been considered in WIMP detection. Among 
them, we consider the isotope targets, 73 Ge and 131 Xe, which have been used at the most sen- 
sitive detectors, CDMS II and XENON10. The recent estimated experimental sensitivities 

n n 

are 0.6 event/kg/day for XENON10 |14j and CDMS II [15(, and the planned experimental 
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TABLE I: Current and planned Dark Matter detectors. XENON10 and CDMS II are the most 
sensitive detectors and SuperCDMS is a future detector. 



Experiment 



Recoil energy range 



Target 



Spin 



Mass (Kg) 



XENON10 
CDMS II 
SuperCDMS 



4.5 ~ 27KeV 
10 ~ lOOKeV 
15 ~ 45KeV 



131 Xe 
73 Ge 



73 



Ge 



5.4Kg 
lOOKg 
100 Kg 



sensitivity is 1 x 10 event/kg/day at SuperCDMS |16fl . The recoil energy ranges are listed 
at Table II for the corresponding detectors. 

Fig. 5 shows the WIMP-nucleon cross sections from Eq.(15) with the experimental limits 
in the WIMP mass range 10 ~ 1000 GeV. The predicted curves and the experimental limit 
contours are for pure proton coupling (a n = 0) and pure neutron coupling (a p = 0). The 
predicted curves show a little different from both isotope targets because of the nuclear spin 
effect. The form factors related to the s pin structure functions are followed by Dimitrov et 



a/.[17| for 73 Ge nucleus and Ressell et al. 18| for 131 Xe nucleus, which are consistent with the 
experimental measurements for MDM in the reasonable precision. The 73 Ge and 131 Xe nuclei 
both contain an unpaired neutron, so the sensitivity to pure-neutron coupling is stronger. 
The predictions are under the current discovery limits for the entire mass range and are 
relatively close to the experimental limits. We should not discard the WIMP discovery 
through the magnetic dipole moment in near future. 



V. CONCLUSION AND DISCUSSION 



A Dirac fermion which was supplied as a singlet within the standard model context may 
mainly interact with the ordinary matters through magnetic dipole moment (MDM). It must 
be more general that the candidate couples to the shifted photon (hypercharge gauge boson) 
than coupling to photon directly near electroweak scale 10 ~ 1000 GeV. Assuming that the 
Dirac fermions account for all the observed dark matter, the MDM constrained by the relic 
density may be as large as 10 -18 ~ 10~ 17 e cm. The WIMP with nonzero MDM interacts 
with targets through MDM of the target nuclei and the cross sections are under the current 
experimental limits. Since the preditions are close to the current experimental limits, the 
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FIG. 5: The predictions for WIMP-nucleon cross sections, (T Pt n, as a function of WIMP mass. The 
light dotted lines are the exclusion limits from the XENON10 experiment and the solid lines from 
CDMS II experiment. Also shown by dash-dotted lines are the experimental sensitivities of the 
planned detector SuperCDMS, respectively. 

detection of the WIMP through MDM should be considered be in near future. 

Besides the direct search of WIMP, there are several experiments to search for dark 
matter indirectly, such like neutrino telescopes, 7-ray telescopes and antimatter telescopes. 
They might give the more serious constraints [19|, but the indirect searches for dark matter 
accompany uncertain informations in astrophysics. This scenario, which the dark matter has 
Dirac fermionic nature and mainly interacts with the ordinary matters through the magnetic 
dipole moment, might not be ruled out. 

We have also studied the possibility of WIMP detection with the method of a photon or 
jet tagging in ILC and LHC. Several tens or hundreds of WIMP events are expected, but 
these signals suffer the huge irreducible background come from the neutrino channel. So the 
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expected signals are swamped by the irreducible background. 
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